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Abstract. We study numerically a two-component A-B spreading model (SMK 
model) for concave and convex radial growth of 2d-geometries. The seed is 
chosen to be an occupied circle line, and growth spreads inside the circle (concave 
geometry) or outside the circle (convex geometry). On the basis of generalised 
diffusion-annihilation equation for domain evolution, we derive the mean field 
relations describing quite well the results of numerical investigations. We conclude 
that the intrinsic universality of the SMK does not depend on the geometry 
and the dependence of criticality versus the curvature observed in numerical 
experiments is only an apparent effect. We discuss the dependence of the apparent 
critical exponent \a upon the spreading geometry and initial conditions. 
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1. Introduction 

The noncquilibrium growth and spreading phenomena are common in nature. The 
examples may include such processes as irreversible adsorption or deposition on the 
surface, liquid invasion in porous media, forest fire, crystal growth, evolution of 
damages in mechanical or electrical systems, epidemic spreading, etc. Presently, 
the different varieties of such phenomena arc extensively studied using the numerical 
models |[ ||, ||, |[ || . The most simple models deal with the ideal species of the same 
size (or type) and the key attention is paid to investigation of interrelation between the 
growth mechanism, surface equilibration efficiency, details of interparticle interactions 
and pattern morphology. 

But the ideal systems of monospecies quite rarely occur in nature and practically, 
certain variations of physical or structural properties are always present in such 
systems. The nonideal polydisperse or multicomponent systems are of a great scientific 
and practical interest |§ ] . Recently, the various growth models with two species or 
phases competition [|9|, |1C, |TT| , P^| , |l3| , and the multicomponent Potts growth 
model [[l5], were investigated. The two-component growth model of Saito and 
Mullcr-Krumbhaar Q generalises the Eden model for the case of two different 
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species, A and B, competition (SMK-model). The authors g| observe the criticality 
of the 2d-model for the case of equal growth rates of two components in a half-space 
planar geometry, i.e. for a case when the growth starts from the linear seed line. This 
criticality results from the rules of species attachment at the growing front, which 
exclude the nucleation of A species in a neighbourhood of domain consisting only 
of B species, and vice versa. The model reveals the competitive growth of separate 
domains and their coarsening. In the limit case of large time or height of the front, h, 
the power-law decay is found for Nab domains, 

Nab ~ h~x, (1) 

where Nab is the number of domains and x=2/3 is the critical exponent for 
meandering interfaces p8| . Later on, Vandewalle and Ausloos Jnj have come to a 
conclusion that the criticality of SMK model is not universal and depends on the 
assumed geometry of two-component propagation. They have shown that if the 
growth starts from the central sites of initial AB configuration, the value of the 
critical exponent \ may differ considerably from that observed for the half-space 
planar geometry (x=2/3). It should be noted that the geometry can effect the fractal 
dimension of DLA-like pattern in a closed cavity, as it was observed in . Recently, 
Batchelor et al. |2(| have shown, that for the usual 2d radial Eden model, the growing 
scaling exponent (3 governing the interface roughening exceeds considerably the value 
of 0=1/3, which is typical for the stochastic growth on a planar substrate (in 1+1 
dimensions), and only approaches this value in the asymptotic limit of the very large 
Eden cluster radius. 

The present work studies kinetics of the two-component spreading on curved 
surfaces in assumption of SMK-model. The evolution of domains, their diffusion and 
annihilation takes place on a two-dimensional square lattice, and spreading is either 
confined inside the circle (concave geometry), or takes place outside the circle (convex 
geometry). For the infinitely large radius of the circle, R— >oo, we revive the half-space 
or planar geometry 0, for other limit case of R—>0 we consider the spreading in a 
free space JTcfl . So this modified model allows us to perform more precise control over 
the influence of geometry on criticality. Using the mean field diffusion-annihilation 
approach and simulations for SMK-model, we calculate the density profiles for wall 
boundaries between A and B domains as a function of distance from the seed wall h for 
seed circles of different radius, R. It will be shown that even in the case when criticality 
remains unchanged, the apparent value of the critical exponent of domain wall spatial 
distribution may vary and it depends on the geometry and initial distribution of A 
and B domains on the seed line. 

The remainder of the paper is organised as follows: Section || presents a brief 
description of the model and simulat ion details. Section || is devoted to our results 



presentation and discussion. Section 3.1 contains theoretical discussions of the issue 
of universality for the competitive growth model with curved seed line. Here, we 
use the mean field approach for description of the domain walls coalescence. In the 
section we present the results of a numerical simulation, performed using the SMK 
algorithm of Q in order to test the theory exploited in section 3.1. Section [| contains 
our conclusions. 
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2. Model and Simulation Details 



We have used the simple two-component SMK growth model |9| with the equal growth 
rates of two components. The simulation takes place on the square lattice. In this 
model, same in Eden model, new particles attach the growing cluster only along its 
perimeter, however in the SMK model the probability of filling any empty site of 
perimeter either by A, or by B species, is in proportion to their concentration at 
the surface of propagation front. The growth starts from the circular seed line and 
develops as the growth front propagation for two different cases: 

• Inside circular closed cavity. This is the case of concave geometry, or inside 
propagation (IP) case. 

• Outside the circular closed cavity. This is the case of convex geometry, or outside 
propagation (OP) case. 

The maximal size of studied systems was 2000 x 2000 grid points and inter-grid 
distance was set as d — 1. For each case, 100-500 samples were averaged. 

The simulation was carried out for different ordered initial configurations, 
including case of a maximal density of domain boundaries with short correlation length 
Xd = 1, where single A and B specie alternate, such as: 
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and, also, in the case of lower densities of domain boundaries, with larger correlation 
length (Xd = 3 for the case displayed), when there is an alternation of larger A and B 
species domains with equal length Xd, such as 
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If the initial number of domain boundaries is equal to N^ B , then the initial linear 
density of the domains walls will be equal to p = N^ B /L, where L = 2irR is the length 
of the seed line. For the case displayed as the seed line we have po = p ma x ~ 1 
and 
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(4) 

1/Xd and 



and for the more common case displayed as the seed line (||) we have po 
2nR/Xd. 

Figure [l] present the typical spreading pattern of IP case with the radius of closed 
cavity R=1000 (a) and the pattern of OP case with radius of circular exclusion R=60 
(b). Here, the thin concentric lines are the time snapshots of the front evolution as 
at the moments of each next 300000 particles attachment. The thick lines show the 
AB inter-domain boundaries. In the OP case the simulation terminates when the 
maximal radius of growth front r max reaches the boundaries of system (r max =1000). 
The initial density of boundaries is chosen to be maximal, po = Pmax ~ 1 such as 
displayed in the seed line (||). Note, that for both IP and OP cases we observe the 
mostly quick decrease of the number of interdomain boundaries Nab in initial periods 
of time, near the surface of circular seed line. 



3. Results and Discussion 

3.1. Mean Field Approach 

The domain evolution can be described in the term of annihilations reaction of domain 
boundary If we associate each domain boundary with fantom particle A, then we can 
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consider the domain coarsening as diffusion-annihilation reaction of type A + A — > 
||. For the linear seed line the time evolution of the density p = Nab/L satisfies a 
mean-field rate equation 

dh~ 9 ' (5) 

giving p ~ ft.~ x , where /i is the height of the front, or equivalently, is the time and a 
is the annihilation rate constant. 

We have x=l/2 for the normal Brownian motion |2]], ^2|, |3), and x=2/3 for a 
two-species Eden Q or SMK Q model. 

Using the simulation for the planar geometry, we have estimated the annihilation 
rate constant as 

a = 2.5 ±0.2, (6) 

where the data used in estimation were averaged over 100 different samples. 

In a case of radial geometry the equation (||) should be modified. For this case 
we have additional source of p rate variations due to the shrinkage or dilation of the 
spreading front for the case of growth spreading inside the circle (concave geometry) 
or outside it (convex geometry), respectively. The domain density p at the distance 
r = i? =F h from the centre of seed circle is equal to 

Nab ~ 
P= 2n(R T hY (7) 

where h is the distance from the seed circle, and, here and hereinafter, the upper 
sign (— ) corresponds to the front spreading inside the seed circle (IP case, concave 
geometry) and the under sign (+) corresponds to the front spreading outside the seed 
circle (OP case, convex geometry). Taking the (|^) into account we obtain instead 
(|J) the following rate equation for the above-mentioned cases of radial spreading with 
p(0)=po as initial condition: 

Introducing the new scaled variables y—pR x , and x—h/R we can rewrite (||) in 
scaled form 

d y i+i/x V tn\ 

— = -ay L+L < x - — — . (9) 

ax 1 =p x 

Equation (||) is a Bernoulli equation and substitution of t^y^ 1 ^ reduces it to 
the differential equation of the following type 

ax 1 =F x x 

e solution of 
form 



The solution of (|To|) with initial condition of t(0) = to = p /R has the following 



t{x) = to(l T x) 1/x ±V(1 T x) 1 - 1 ^ + 1 T V , (11) 

where 

V=a/(t (l - x)) = aRpl /x /{\ - x )- (12) 
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Taking into account the inverse substitution y=t x , we obtain the following final 
solution of differential equation (||) 

y ^ ~ (i t x) [±v(i i x y- l /x + 1 T v] x ' (13) 

where yo=poR x . 

For the IP case, when the spreading takes place inside the circle, the value y{x) 
goes through the minimum at the point 

'-(t^TTTjV" 1 for V> 1/(1 -x). 



Vi-V^y - /v (14) 

for V< 1/(1 -x). 

For the case when V — ► oo, we have 

s m in = hm x min (y)=l - X^ 1 , (15) 

and 

y m in= lim y(x m in)=a x x 1 ~ x - (16) 

Figure || presents the y/y m in versus x dependencies obtained from (|l3|) , (|l6|) for 
X=2/3 at different values of V for the case of growth spreading inside the circle (solid 
lines) and outside the circle (dashed lines). The slope value -2/3 accounts for the slope 
of lines in the limit V — > oo for the spreading both inside and outside the circle, and 
the slope value -1 accounts for the slope of the lines in the limit of x — > oo for the 
growth spreading outside the circle. 

From ([!]), we can define the apparent value of the critical exponent Xa as 

d\nN AB 

Xa = -^h— (17) 
Figure ^ shows a plot of the apparent critical exponent Xa versus scaled distance 
from the seed surface x for both concave (dash lines) and convex (solid lines) 
geometries. We see that for the propagation started from the curved seed line the 
value of the apparent exponent Xa is not constant and strongly depends on x. 

Not far from the seed surface, at x — > and high initial density po = p ma x ~ 1 
and in the limit V — > oo, the intrinsic value Xa — X — 2/3 is observed. The observed 
Xa values substantially decrease with initial density decrease. With increase of x, the 
value of Xa always increases for the IP case. For the OP case, the Xa goes through 
the maximum. We should stress that these are only apparent changes of critical 
exponents, estimated on the basis of (0), because all calculations were done for the 
universal system in assumption of constancy of intrinsic value (x = 2/3). 



3.2. Numerical Simulation 

The dependencies of y=pR x upon x=h/R obtained as results of numerical simulations 
are depicted in figure |] for the IP case (a) and OP case (b). The different points 
correspond to the different radii of initial circular seed lines. All these data are 
obtained for the case of the maximal initial density of domain boundary po — p max ~ 1 
with ordered configuration of A and B the species, as in the seed line (||). The 
continuum limit for this case corresponds to the situation when R 3> 1 and we can 
neglect the lattice discreteness. The solid lines in figure || obtained from (|l3|) for 
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the case of x=2/3, V — > oo in the continuum limit. We see that in this limit the 
coincidence between the results of the mean field approximation and the computer 
simulation is rather good. It is important to note that in the limit of higher initial 
density of AB boundaries, all data obtained at different values of R in the scaled 
co-ordinates pR x versus x=h/R fall on the same curves for both IP and OP cases. 

For the relatively small values of initial density po, when we start propagation 
from configuration with limited number of domain boundary Nab ~ 10 — 100, the 
simulation results can deviate substantially from predictions of the mean-field theory 
as it is shown on figure || for the IP case (a) and for the OP case (b) . The reason of 
such deviation is unclear by far but we think that it simply reflects the limitations of 
the mean field approach, which follow from neglecting by the long range correlations 
between the AB domain boundaries. 

4. Conclusions 

We may conclude that the apparent critical exponent \a which is determined through 
treatment of the numerical experiments results for SMK model on the basis of (|l7|), 
is not universal or constant value. For growth spreading in the case of nonplanar 
geometry the apparent critical exponent \a depends on the height of the AB front, 
initial numbers of AB domain and type of the spreading geometry. This conclusion 
remains true for SMK model even in the case when the model preserves its internal 
universality as defined by the classical critical exponent \ = 2/3. The most good 
agreement between results of numerical simulation and of the mean field approximation 
is observed only at high initial number of boundaries N^ g . At low initial number of 
boundaries N^g, the mean field approximation allows to reach only the qualitatively 
true description. In this case, the absence of quantitative description may be explained 
by the effect of long-range correlations, which result in the faults of the mean field 
description. 

We believe, that results obtained in the present work coincide qualitatively quite 
well with data presented in |l(J. The authors of this work have shown that for 
initial configuration consisting only from two species AB the number of boundaries 
is approximately conserved at the level of Nab ~ 2 and in this case they obtained 
Xa ~ 0. Assuming for propagation under such conditions N^ B = 2,-R~l,X = 2/3 
and the value a = 2.5 ± 0.2 (g) we have V « 1 as estimated from ([l2|) . As it follows 
from our data presented at figure ||, when the V values are so low, the \a value is 
small and practically equal to zero at high enough distances from the seed centre, 
which is in correlation with data obtained in p0| . We believe that such apparent 
alterations of criticality or scaling exponential functions may occur in other models 
of radial growth, as it was observed, e.g., in |^(| for Eden model. The similar tasks 
are of interest for further investigations with the extended range of models, number 
of system components, space dimensionality, etc. 
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Figure 1. The typical growth spreading inside the circular cavity of radius 
R=1000 (a) and outside the circular exclusion of radius R=60 (b) for two-specics 
SMK model. The thin conccntrical lines are the time snapshots of fronts for each 
next group of 300000 particles attachment. The thick lines show the AB domain 
interfaces. 
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Figure 2. The mean field dependence of the scaled density y/y m in {y = pR x ) 
versus scaled height of the front x = h/R for the growth spreading inside the 
circular cavity (solid lines) and outside the circular exclusion (dashed lines) at 
different values of V=aRpQ — x)- The value of y m i n as defined by (|l|) 

corresponds t o th e minimum of y(x) function in the point x = x m i n = 5/9 (for 
X = 2/3, see (|15|)) for growth spreading inside the circular cavity in the limit of 
V -> oo. 




Figure 3. The effective critical exponent \ a versus scaled distance from the seed 
surface x = h/R for both concave (dash lines) and convex (solid lines) geometries 
at different values of V=aRp^ x /(l — Arrow shows the classical value of 
X = 2/3 for the planar geometry. 
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Figure 4. The scaled density y = pR x versus scaled height of the front x = h/R 
for the growth spreading inside the circular cavity (a) and outside the circular 
exclusion (b). The different points are the simulation results for different cavity 
(a) or exclusion (b) radiuses, which are presented on the figures. All data are 
obtained for the limit case of po = Pmax, for configuration shown by line (Eh. 
The solid lines are obtained from the mean field equation (jllj) where we used 
a = 2.5 ± 0.2 flfj) estimated by simulation for the case of planar geometry. Point 
of minimum x — x m i n — 5/9 « 0.561 is defined by the (M) for x = 2/3. 
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Figure 5. The scaled density y = pR x versus scaled height of the front x = h/R 
for the growth spreading inside the circular cavity or radius R = 800 (a) and 
outside the circular exclusion or radius R = 60 (b). The solid irregular lines are the 
simulation results for different initial numbers of domain interfaces (a): N^ B = 
20(1), 40(2), 80(3), N ma *(4), and (b): N% = 10(1), 20(2), 40(3), 90(4), N m ax(5)- 
Dashed lines are obtained for the same numbers of domain interfaces through the 
mean field result ( ]l3[ ) where we used a = 2.5 ± 0.2 (H), estimated by simulation 
for the case of planar geometry. 



